An Upper Bound for the Spin- Wave Spectrum 
of the Heisenberg Antiferromagnet|^ 



> 

i> 
o 
\o 
o 

B 



Tsutomu Momoi 



Department of Physics, University of Tokyo, Hongo 7-3-1, Bunkyo-ku, Tokyo 113, Japan 



Abstract 

We study the spin-wave spectrum of the Heisenberg antiferromagnet on a 
bipartite lattice. The spin-wave spectrum on a Neel-ordered ground state is 
bounded as 6{k) < c\k\, where c gives an upper bound for the spin-wave velocity. 
In the large- S" limit, the upper bound c coincides with the result of the spin- wave 
theory. 
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. Ground-state properties of the quantum Heisenberg antiferromagnet have been stud- 

ied extensively. The existence of Neel order in the ground state has been proven for the 
three-dimensional spin-S* Heisenberg antiferromagnets on bipartite lattices.0 " & In the two- 
dimensional systems, the occurrence of symmetry breaking at T = has been proven for 
higher spin (S > 1) models and for anisotropic models.i"& In the spin 1/2 Heisenberg 
antiferromagnet, many theoretical and numerical studies indicate that the system has a 
long-range order, though there is no rigorous proof (to date). 

The Nambu-Goldstone theorem for quantum spin systemsi-* states that, if there is a 
Neel order, there are gapless elementary excitations. In the Heisenberg antiferromagnet, 
the spin-wave theory!' 0^ succeeded in giving a precise description of the Goldstone bosons, 
^To be published in J. Phys. Soc. Jpn. 63 (1994) No. 7. 
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i.e., magnons. The linear spin-wave theory predicts the gapless /c-hnear dispersion relation, 
s(k) ~ 2VdSJ\k\, for small k. 

In this letter we consider the quantum Heisenberg and XX Z antiferromagnets on the 
(^-dimensional Lx ■ ■ ■ x L hypercubic lattice A C Z*^ and we consider the case oi d > 2. The 
Hamiltonian is given by 

Hk^J (S^S^ + S!S] + AS^S^), (1) 

where < A < 1 and the summation runs over all the nearest-neighbor sites. The system 
size of A is A^(= L'^). We show an upper bound for the spectrum of elementary excitations 
on the Neel-ordered ground states. We find that the excitation spectrum is bounded as 
e{k) < c\k\, where the constant c gives an upper bound for the spin-wave velocity. For 
S ^ 1, the upper bound c coincides with the result of the spin- wave approximation. In the 
Appendix, we show some rigorous results for the large-S" Heisenberg antiferromagnet. 

Let us consider the system under the small staggered magnetic field whose Hamiltonian 

is 

HAiB) ^Ha- BMa, (2) 

where 

MA = T.Sf-i:S!. (3) 

ieA ieB 

We denote the normalized ground state as |$gs,b)- We take a thermodynamic limit applying 
an infinitesimally small field B. The staggered magnetization is given by 

ms = lim lim ^(^gs.bIA^aI^gs,^)- (4) 
As an excited state, we consider the following standard trial state 

\^B{k)) ^ S^,\^GS,B)m\^GS,B)\\, (5) 

where — N~^/'^ J2i S- exp{ik ■ ri), k — {ki,k2, ■ ■ ■ ,kd) and 

||'5'||$gs,b)|| = {^GS,B\^-k^k\^GS,By^'^- When spins lie in the xy plane, the operation of Sf 
flips the spin at the site i. The excitation energy of is given by 

=limlim{(*5(fc)|7iA(S)|*5(fe)) (6) 

BIO Ajoo 

-($gs,b|^a(5)|$gs,s)}. 
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We take the momentum as fc 7^ 0, 7^ (vr, . . . , vr) and = 2711^/ L (0 < < L — 1), 
so that |$gs,b) and B{k)) are orthogonal This trial state is called the Bijl-Feynman 
single-mode approximation.0' A similar trial state was studied by Horsch and von der 
Linden in the two-dimensional system with no magnetic field.El^ Using the state @, we 
show an upper bound of the spin-wave spectrum as follows.^ 



Theorem: If the ground state has a Neel order, i.e., if mg > 0, the energy spectrum 
e{k) is bounded as 

e{k) < + Py) ^^^(1 + A^^) + ^^(^ + ^,)yr^, (7) 

where 

"""-ISfe]^**--! ^ SfS?l*GS,B) (8) 

(«j>eA 



for a = X, y and z, and 



1 



^ i=i 



Proof. As in ref. 0, the excitation energy of \'^B{k)) is calculated as 

^\1'GS,B|'J-fe'Jfc|1'GS,B/ 

2J(1-7,)(<I>gs,b| E i-S!S^ - SfS])\^GS,B) + B{^gs,b\Ma\^gs,b) 

('•,i>GA 



2iV(<l>GS,B|5ife5||<l>GS,B) 

The excitation energy in the thermodynamic limit is given by 



(10) 



where >S'j_(fc) denotes the structure factor of the Neel-ordered ground state, 

Slik) = limlim($GS,B|^!.^fc1$GS,B), (12) 

BIO Ajoo 



^After the submission of this manuscript, we learned that Stringari had previously obtained a better 
upper bound of the spin-wave spectrum.© We thank Mr. K. Totsuka for letting us know Stringari's paper 
(ref. ^) . New points of the present paper are that we obtain both upper and lower bounds of the excitation 
spectrum of the Feynman state and that we give rigorous arguments on the large- limit. 
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which describes correlations of spins transverse to the magnetic field. To bound the structure 
factor for the Heisenberg model from below, we use the correlation inequality by Sha=try,0) 



2Sl{k) > ^ 

VPx + PzVl + Ik 

This bound can be extended to the anisotropic model (H) in the form 

msVl - Ik 



(13) 



2Sl{k) > 



Using eqs. ([TT|) and (0), we obtain eq. 



p^(l + A7fc) +p,(A + 7fc) 



(14) 
■ 



Consequently, the lowest spin-wave spectrum is bounded from above by a gapless fc-linear 
dispersion relation. We thus find that the spin-wave velocity Vs is bounded from above in 
the form 

< V^2c/(1 + A) (p^. + p,)(p^ + Py)J/mJ^. (15) 

The expectation value e{k) can be bounded from below as well. The transverse structure 
factor S^{k) is bounded from above in the form0 I-* 



2Sl{k) < 



{Px + Py)il ~ Ik) 



1/2 



(16) 



A(l+7fc) 

Note that theorem 4.1 of ref. |l| is applicable to the Hamiltonian (2). Using this inequality, 
we obtain 

' (17) 



6{k) > 2rfJ^A(p, + pJ(l-7fc2). 

(We remark that this gives only a lower bound for the expectation value of the single-mode 
approximation and not for the spin- wave spectrum.) Both the upper and lower bounds have 
the /c-linear dispersion relation for small k. We thus find that e{k) has the gapless fc-linear 
dispersion relation. 

Finally, we discuss the large- S* limit of the spectrum. For the Heisenberg antiferromagnet, 
we have 



lim rUs/ S 

S^oo 

lim p^/S"^ 

b~*oo 

lim py/S^ 



lim pJS^ = 0. 

S— >oo 



(18) 
(19) 
(20) 
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We show a proof in the Appendix. From these rigorous results, we find that in the large-S* 
hmit the upper bound for the Heisenberg model coincides with the lower bound and we 
obtain 

eik) = 2dSJ^l-^k^. (21) 

This is the same result as that given by the spin- wave theory.!' 0* We thus find that the 
single-mode approximation gives a precise spin- wave spectrum in the large- limit. For 
the large-S* XX Z model (0 < A < 1), it is expected that the ground state becomes 
the Neel state and we have px — S^, py ^ 0{S), pz — 0{S) and nis ^ S, though we 
have not proven it. Then we expect that the upper bound for the spectrum behaves as 
2dJS\J (1 + A7fe)(l — 7fe). This spectrum coincides with that given by the linear spin-wave 
theory for the XX Z model.0^ 

The author would like to thank Professor M. Suzuki for helpful discussions and for 
critically reading this manuscript. The author is also grateful to Dr. T. Koma and Profes- 
sor H. Tasaki for useful comments. 



Appendix. Some rigorous results for the large-*? Heisenberg antiferromagnet 

In this appendix we show bounds for nis, Px, Py and pz in the Heisenberg antiferromagnet 
on square and cubic lattices. From these bounds, we obtain the S dependence of these values 
in the S ^ oo limit. 

The quantities rris, Px, py and pz are bounded from above by the norms of local operators. 
Hence, we obtain 

m^<S (22) 

and 

Pa<S'^ {a = x,y,z). (23) 

The value {—px — Py — Pz) is the ground-state energy per bond, whose bounds were shown 
by Anderson0) in the form 

< {p. + Py + Pz) < S{S + l/2d). (24) 
To bound from below, we use Koma and Tasaki's inequalityB* 

ms > 3 lim — ($Gs,iJ=o|A<A^|$GS,i3=o) (25) 



and the inequality by Dyson, Lieb and Simon& and Jordao Nevez and Fernando Perezi^ 



lim ^($gs,b=o|-Ma^|'^'gs,b=o) 

> SiS + l)/3 - (p, + + p,y/^h/V6, 



(26) 



where 



1/2 



{2nY\l+jkJ ■ ^^^^ 

The integral Id is finite for d > 2. It was shown that the r.h.s. of the inequality (p6|) is 

positive for S* > 1 in the two- and three-dimensional systems. Using eqs. (p^)-(l26|), we 

obtain bounds of in the form 



S{S + 1) - ^3S{S + l/2d)Id/V2 

< < S 



1/2 



(28) 



for 5* > 1. Thus we find 



lim rris/S = 1. 

5— >oo 



(29) 



To bound px, Py and pz from below, we use eqs. (|T^) and (|T^), from which we find 



From symmetry of the Hamiltonian, i.e., Py = Pz, we obtain 



(30) 



rris <Px + Py 



(31) 



Using eqs. ( PTf ) and (P^), we obtain 



ms' + Pj, < 5(5+l/2rf), 



(32) 



where we have used Py = pz- Using eqs. (p3|), (^) and (p2|), we obtain bounds for po, 
(a = X, ?/, z) in the form 



2m^^ - S{S + l/2(i) < - Py < Px < 5l 



(33) 



Using eqs. (|33D and (p9|), we find 



lim p^/S"^ = 1 



(34) 



and 



lim Py/S"^ = lim Pz/S'^ = 0. 

S^co S^oo 



(35) 
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